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ABSTRACT 

We report on partial and global generalized synchronization in exemplar bi-directionally coupled 
non-identical chaotic systems with multiple time delays. We derive conditions for such synchro- 
nization regimes. The general approach is applied to the case of multiple time delay laser systems 
with electro-optical feedback. We also study the effect of signal sampling on the signatures of time 
delays in the autocorrelation function of the multiple time delay laser output. 

1. INTRODUCTION 

The synchronization of systems of coupled oscillators is important in many disciplines of science. 
In a chaos based secure communication scheme synchronization is vital for message decoding [1]. 

Recently, delay differential equations (DDE) have attracted much attention in the field of 
nonlinear dynamics. The high complexity of the multiple time-delayed systems can provide a 
new architecture for enhancing message security in chaos based encryption systems [2]. In such 
communication systems message decoding would require chaos synchronization between multiple 
time-delayed transmitter and receiver systems. 

Several kind of synchronization can arise in interacting systems. Generalized synchronization, 
first studied for uni-directionally coupled systems, is defined as the presence of some functional 
relation between the states of response and drive, i.e. y(t) = F{x(t)) [3]. As underlined in [4], for 
secure chaos based communication schemes the use of generalized synchronization mode could be 
preferable in comparison with the non-generalized modes of synchronization. Recently, generalized 
synchronization between bi-directionally linearly coupled non-identical Lorenz and Rossler mod- 
els has been investigated numerically in [5] (by Zhigang et at), where the concept of generalized 
synchronization was extended to mutually coupled systems. It was shown that with increasing 
coupling strengths the system changes from partial to global generalized synchronization. 
It is noted that bi-directionally coupled time delay systems with a single time delay were stud- 
ied previously in by Heil et al [5]. Paper by T. Heil et al presents experimental and numerical 
investigations of the dynamics of two device-identical optically coupled semiconductor lasers ex- 
hibiting a delay in the coupling. They find well-defined time lag between the dynamics of the 
two lasers and an asymmetric physical role of the subsystems. Further, in paper by Fischer et al 
[5] the authors show that isochronous synchronization between delay coupled oscillators can be 
achieved by relaying the dynamics via a third mediating element. In other words they consider 
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three bi-directionally coupled oscillator in a line configuration and demonstrate that zero-lag syn- 
chronization can occur over long distances through relaying. 

To the best of our knowledge partial and global generalized synchronization in time delayed sys- 
tems has yet to be investigated. 

This paper is devoted to the study of partial and global generalized synchronization between 
mutually coupled chaotic systems with multiple time-delays. Wc find analytically the stability 
condition of the global generalized synchronization mode between bi-directionally coupled multi- 
ple time-delayed systems and apply the approach to the multiple time delay semiconductor lasers 
with electro-optical feedback. 

In a time-delayed-chaos-based communication, identification of time delays can allow the eaves- 
dropper to extract the message successfully using a simple local reconstruction of the time delay 
system. The delay time can be revealed by a number of means including via the autocorrelation 
coefficient [6]. With this in mind, in this paper we also investigate the effect of sampling rate 
on the autocorrelation coefficient of the laser system with electro-optical feedback. The principal 
result of this study is the demonstration that the recovery of the delay times from the autocorre- 
lation coefficients of the laser output depends on the sampling rate of the signal. 



II.GENERAL APPROACH 

We consider synchronization between the bi-directionally coupled double-feedback systems of gen- 
eral form , 

dx 

— = -aix + mifixri) + m^fixr^) + Ki{y - x), (1) 

dy 

— = -a2y + msfivrs) + mj{yrj + K2{x - y), (2) 

where / is differentiable generic nonlinear function; q;i,2 are the relaxation coefficients for dynam- 
ical variables x and y; mi, 2 and m3,4 are the feedback rates for systems x and y, respectively; ti,2 
and r3 4 are the feedback delay times in the coupled systems; Ki^2 are the coupling rates between 
X and y. Throughout this paper Xr = x{t — r). 
We find that under the conditions 

cci = a2, mi = 777,3, "72 = 1714, K2^ K, ti = T3, T2 = T4 (3) 

Eqs.(l) and (2) admit the synchronization manifold 

y = x. (4) 

This follows from the dynamics of the error A — x — y (with conditions (3)) 

^ = -«iA + miArJ'ixr,) + m2ArJ'{xr,) - 2KA. (5) 
dt 2 



obtained by use of a Taylor series. Here f'{x) denotes for derivative of / and the derivative should 
be bounded. Solution y = x under the conditions (3) is also follows from the symmetry of Eqs.(l) 
and (2). According to the Lyapunov-Razumikhin functional approach (see, e.g. [2] and references 
there-in), the sufficient stability condition for the synchronization regime y = x {4) as: 

ai + 2K > |mi(sup/'(a;^J)| + |m2(sup/'(a;^))|. (6) 

Here sup f'{x) stands for the supremum of /' with respect to the appropriate norm. 

Next we consider the effect of parameter mismatches on synchronization . For the parameter 
mismatches, ai ^ a2, mi ^ m^, m2 7^ ^4, Ki 7^ A'2, ti 7^ rs, T2 7^ r4 is clear that complete synchro- 
nization, Eq. (4) is no longer the synchronization regime for coupled systems (1) and (2). Then for 
such a case one can try using the auxiliary system method to detect generalized synchronization 
between bi-directionally coupled x and y [5]. According to [5] bi-directionally coupled systems x 
and y can be considered as two uni-directionally coupled systems driven by x and y. Then for each 
of the uni-dircctionally coupled systems one can apply the concept of generalized synchronization: 
that is given another identical driven auxiliary system z{t), generalized synchronization between 
x(t) and y{t) is established with the achievement of complete synchronization between y{t) and 
z{t). Thus, in the case of bi-directionally coupled systems for generalized synchronization one 
has to investigate the possibility of two complete synchronization manifolds, one for x driving, 
and another one for y driving. Depending on the achievement of complete synchronization due to 
the coupling strengths between the systems, partial (when the complete synchronization occurs 
between the systems driven either by x or y) and global (when complete synchronization occur 
between the systems driven by both x and y) generalized synchronization can arise. Applying the 
error dynamics approach to investigate the complete synchronization between the systems (driven 

by y) 

dx 

— = -aix + mif{xr^) + 1712 f{xr^) + Ki{y-x), (7) 

— = -aixi + mifixi^n) + m2f{xi^r2) + Ki{y - xi), (8) 
and the systems (driven by x) 

dy 

— = -a2y + msfiyrs) + nnfiyrj + K2{x - y), (9) 



Jl = -a2yi + msfiyi^rs) + m4f{yi,rj + K2{x - yi), (10) 



dyi 
dt 

we find the sufficient stability conditions for the complete synchronization x — xi 

Ki > |mi(sup/'(a;ri))| + |m2(sup /'(x^-J)! - ai. (11) 



3 



and for y — yi 

K2 > |m3(sup/'(y^))| + |m4(sup/'(y^J)| - as- (12) 

Then for (i^Ti, A'2) > max(7^i,7^2) complete synchronization occurs, which means that we have a 
case of global generalized synchronization between the bi-directionally coupled systems x and y. 
We notice that the stability conditions (6), (11) and (12) derived using the Lyapunov-Razumikhin 
approach is a sufficient one: it assures a high quality synchronization for a coupling strength 
estimated from the stability condition, but does not forbid the possibility of synchronization with 
smaller coupling strengths. The threshold coupling strength can be estimated by the dependence 
of the maximal conditional Lyapunov exponent A of the error dynamics on coupling strength 
K:i.e. from \{K) =0 [7]. Suppose that and K2 are the value of coupling strengths found from 
\{K) ~ for auxiliary systems (7-8) and (9-10), respectively. Then for {Ki,K2) < min(ii'", ) 
there is no synchronized state, for mm{K'^ , K2) < {Ki,K2) < max(irf , ^2 ) there is partial and 
for {Ki, K2) > max(ir", K2) there is global generahzed synchronization. 



III.WAVELENGTH CHAOS MODEL 

Synchronization of chaotic semiconductor lasers is of significant practical importance, as these 
lasers have potential application in high-speed secure communications [1]. The laser system con- 
sidered in this paper is an electrically tunable Distributed Bragg Reflector (DBR) laser diode with 
electro-optical feedback. This system was proposed in [8] as a chaotic wavelength signal generator 
for chaos based secure communication. The wavelength of the chaotic carrier is described by the 
following dynamical equation: 



^ "^^^^ + si"'(^A(t - r) - $0) (13) 



where A is the wavelength deviation from the center wavelength Aq;D is the optical path difference 
of the birefringent plate that constitutes the nonlinearity;<l>o is the feedback phase;r-thc feedback 
loop delay time; T is the time response in the feedback loop;/?^ is the feedback strength. With 
X — and m — ^^i^ we rewrite Eq.(13) in the following normalized form: 



^ -ax{t) +msin'^{xr - $o)- (14) 



dt 



Notice that in (14) we have scaled the time with aT{a is the relaxation coefficient). In the fol- 
lowing we consider chaos synchronization between Bi-directionally linearly coupled electro-optical 
semiconductor lasers with double feedbacks. 



IV.BIDIRECTIONALLY COUPLED SYSTEMS 
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The lasers to be synchronized are described by the following equations 

ax{t) + mi sin^(xT-i - $o) + '^2 sm'^{xr2 - ^o) + Ki{y - x), (15) 

ay{t) + rris sin^(y^3 - $o) + m4 sin^(y^4 - $o) + K2{x - y), (16) 

Applying the general approach from Section II we find that systems (15) and (16) can be syn- 
chronized on the synchronization manifold (4) under the existence conditions (3). The sufficient 
stability condition for the synchronization manifold y — x ior the coupled systems (15-16) is : 

2K > (mi + 1712) - a. (17) 

Considering the effect of parameter mismatches on the synchronization regime by use of the 
results of Section 11, we establish that the global generalized synchronization between x and y 
occurs if the coupling strengths {Ki, K2) > meix{Ki, K2) where Ki and K2 satisfy the inequalities: 

Ki > {mi + m2) - a, (18) 

> (ma 7714) - a. (19) 

Under the condition (18) complete synchronization occurs between the y driven auxiliary systems 
and (19) is the sufficient stability condition for the complete synchronization between the auxiliary 
systems driven by x. 



dx{t) 
dt 

dyjt) 
dt 



V.NUMERICAL SIMULATIONS AND DISCUSSIONS 

Numerical simulations fully support the analytical results. In the simulations significant attention 
is given to the study of partial and global generalized synchronizations between non-identical sys- 
tems. This is due to the possible advantages of generalized synchronizations in secure chaos based 
communication schemes, as underlined in [4]. With the security issue in mind, we also investigate 
the effect of the sampling rate on the autocorrelation function modulations. 
Numerical modelling of Eqs.(15) and (16) were conducted using DDE23 program in Matlab 
(R2008b). Below we present results of numerical simulations for symmetrical couphng {K — 
Ki — K2) to emphasize the fact that, despite the symmetry in the coupling in bi-directionally 
coupled systems there is asymmetry in achieving synchronization between these directions. 
We have also calculated the cross-correlation coefficients [9] using the formula 

C(A«) = <{.(t)-<x»(y(t + At)-<y»> ^^^^ 

^< ix{t)- < X >)2 >< {y{t+At)- <y>y> 
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where x and y are the outputs of the lasers, respectively; the brackets< . > represent the time 
average; At is a time shift between laser outputs. This coefficient indicates the quality of syn- 
chronization: C=l means perfect synchronization. 

Figure 1 shows complete synchronization between mutually coupled identical laser systems 
with double time delays. 

Before considering synchronization between non-identical systems, we would like to dwell on the 
one aspect of the security for chaos based communication schemes, namely the effect of the sig- 
nal's sampling rate on the autocorrelation function. Figure 2 demonstrates the autocorrelation 
coefficient (here the autocorrelation coefficient will be denoted by Ca and obtained from equation 
(20) when x — y ) ior the output of x laser, Eqs.(l-2) for double time delays, with n = 3, T2 = 5, 
mi = 2.3, m2 = 2.5, $0 — 7r/4, a = 1. It is clearly seen that time delays can be easily recovered 
from the autocorrelation coefficient, as it exhibits extrema at values of the time delays or their 
multiples and combinations. In this connection we would like to emphasize the following important 
point. The number of data points i.e. sampling rate used to calculate autocorrelation coefficient 
is highly significant to detect the modulations of the autocorrelation coefficient at time delays. 
In other words the sampling rate is of importance for the security of chaos based communication 
systems. Indeed for figure 2 the number of data points used was 65, which allowed recovery of the 
time delays. In figure 3 the number of data points was 25. Figure 4 contains 5000 data points. It 
is noted that both undersampling (fig.3) and oversampling (fig.4) do not allow extraction of the 
correct time delays from the modulations of the autocorrelation coefficients. 

In the remainder of the paper we focus on the synchronization between non-identical sys- 
tems with double time delays. In figures 5-7 we present the results of numerical simulations of 
Eqs. (15-16) for for $0 = 7r/4,mi = 20, 1712 = 15, = 4, 777,4 = 3, ri = 2, T2 = 1, Ts = 6, T4, = 7. Fig- 
ure 5 demonstrates the unsynchronized state between non-identical lasers for coupling strengths 
Ki — K2 — 0.1. It is noted that for such couphng strengths the error signal between the auxihary 
systems in both directions does not approach zero. In figure 6 partial generalized synchronization 
between the interacting lasers is presented for the coupling strengths Ki = K2 = 3.2. In the case 
of partial generalized synchronization the error signal in one direction (7/1 — y) approaches zero 
with time, i.e. synchronization is achieved between the auxiliary systems driven by laser x. In the 
reverse direction there is no synchronization,i.e. the synchronization error signal Xi — x does not 
tend to zero with time. 

Figure 7 presents the case of global generalized synchronization between mutually coupled 
lasers x and y for Ki = K2 = 12. In this case the auxiliary systems driven by both lasers are 
synchronized, i.e. both error signals yi — y and Xi — x tend to zero. We notice that in the case of 
partial generalized synchronization, despite the symmetry in the coupling, the coupling directions 
are not equivalent; there is synchronization in one direction and no synchronization in the reverse 
direction. In other words there is always a master laser and a slave laser, despite the symmetrical 
coupling. With increasing coupling strengths this asymmetry in the status of lasers is eliminated: 
both systems x and y control each other in the sense that both driven systems (by x and y) are 
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synchronized and eventually x and y enter a state of global generalized synchronization. 
It is noted that the range of coupling strengths for the asymmetry is determined by the values 
of K when the maximum conditional Lyapunov exponents for the auxiliary systems reach zero. 
Figure 8 depicts the maximum transversal (conditional) Lyapunov exponents A!^(A) and (<)) 
versus the coupling strength for the auxiliary systems Xi and x driven by y and yi and y systems 
driven by x, respectively. As evidenced from the numerical simulations, for the data used in our 
simulations the value of the maximum transversal (conditional) Lyapunov exponent reaches zero 
A?" ~ at ~ 3.1 (when the auxiliary systems yi and y driven by x laser is synchronized) and 
A^ ~ at K2 ~ 9.2 (when the auxiliary systems Xi and x driven by y laser are synchronized). 
As suggested in [5] the more chaotic system plays the role of a master, and less chaotic systems 
assume the role of the slave. The simulation results thus support the fact that x laser is more 
chaotic, and hence it behaves as the master laser, as the largest Lyapunov exponent for x laser is 
^xmax ~ 0.93, and for y laser Xymax ~ 0.10. It is noted that x laser is more chaotic than y laser 
because of higher feedback values. 

Finally we consider the role of parameter mismatches to achieve synchronization between 
uni-directionally linearly coupled non-identical time-delayed systems and show that parameter 
mismatches are of crucial importance to achieve synchronization. 



Let there be certain mismatches between the relaxation coefficients:ai = a — S, and a2 = ct + S. 
Then by investigating the dynamics of error — ?/ = A it is straightforward to establish that 
synchronization y = Xr^ is possible under existence conditions mi = 777.3,777.2 = 772.4, q;2 — oti = K 
and is stable if 0:2 > l^sl + |^4|- 

Detailed study estabhshes that independent of the relation between the delay times in the coupled 
systems and the coupling delay time, only retarded synchronization with the coupling delay time 
is obtained (figure 9). Most importantly with or without parameter mismatch neither complete 
nor anticipating synchronization occurs. 

We mention that, for example in the case of non-linear (sinusoidal) coupling for identical drive and 
response systems with electro-optical feedbacks, depending on the relation between the feedback 
delays and the coupling delay time retarded, complete or anticipating synchronization can occur 
[10-11]. These results are of significant interest in the context of relationship between parame- 
ter mismatches, coupling forms and synchronization. Indeed, having in mind possible practical 
applications of anticipating chaos synchronization in secure communications, in the control of 
delay-induced instabilities in a wide range of non-linear systems, etc, by choosing the " appropri- 
ate" parameter mismatches and coupling forms certain types of synchronization can be achieved. 




(21) 



dy{t) 



a2y{t) + 7773 sin^(7/^3 - $0) + "^4 sin^(7/^4 - $0) + Kx. 



(22) 



dt 
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VI.CONCLUSIONS 



We have investigated partial and global generalized synchronization regimes between two bi- 
directionally coupled time-delayed systems with multiple delays. We have derived the conditions 

for such synchronization regimes. It is well known that in chaos based communications chaos 
synchronization between the transmitter and receiver is vital for the message decoding. As such 
synchronization conditions are of importance for communication between the systems. We have 
tested in the case of laser systems with electro-optical feedback. The results are important from 
the point of view of transition from the disordered state of interacting complex systems to the 
partial and global order states. These findings are also of importance for secure chaos-based 
communication systems in the context of advantages of generalized synchronization mode in com- 
parison with the non-generalized regimes of synchronization. 

Also by studying the effect of the sampling rate on the signatures of time delays in the autocor- 
relation function of the multiple time delay laser output we have demonstrated that the optimal 
sampling rate is essential to recover the time delays. This property is also important for secure 
chaos-based communication systems. 

Acknowledgements.-This research was supported by a Marie Curie Action within the 6*'* Euro- 
pean Community Framework Programme Contract. 
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Figure captions 



FIG. 0. Schematic experimental set-up for wavelength chaos synchronization between the trans- 
mitter and receiver laser diodes with electooptical feedback: DBR LDl and DBR LD2 are the 
Distributed Bragg Reflector transmitter and receiver laser diodes, respectively; BS, beamsplit- 
ter; M, mirror; BP, birefringent plate between crossed polarizers(not shown in the figure); PD, 
photodiode; DL, delay line; LPF, low-pass filter; G, optoelectronic gain; 1, DBR-section injection 
current; 01, optical isolator to provide unidirectional coupling between laser diodes. 
FIG. 1. Numerical simulation of mutually coupled identical electro-optical laser systems, Eqs.(15- 
16). Complete synchronization: time series of x{t) (solid line) and y(i) (dashed line) for a = 
l,$o = vr/4,mi = = 5,m2 = 1714 = 7, Ki = K2 = 9,ri — ts — 2,T2 — T4 — 6;G is the 
correlation coefficient between x and y. Dimensionless units. 

FIG. 2. The autocorrelation coefficient Ca of a; laser output for parameters, Eqs.(15)(i^ri = 0) for 
Ti — 3,T2 — 5, mi = 5, 7712 — 7, a — 2,^0 — tt/A. The number of data points N=:65; Sampling 
time=300/N. Dimensionless units. 

FIG. 3. The autocorrelation coefficient Ca of x laser output for parameters, Eqs.(15)(A'i = 0) for 
Ti = 3,r2 = 5, mi = 5, m2 = 7, a = 2, $0 = 7r/4. The number of data points N=25; Sampling 
time=300/N. Dimensionless units. 

FIG.4. The autocorrelation coefficient Ca of x laser output for parameters, Eqs.(15)(i^i = 0) for 
Ti = 3, T2 = 5, mi — 5,1712 — 7, a — 2, — 7r/4. The number of data points N=5000; Sampling 
time=300/N. Dimensionless units. 

FIG. 5. Numerical simulation of mutually coupled non-identical electro-optical laser systems, 
Eqs.(15-16): Unsynchronized state between x and y for $0 = 7r/4,mi = 20, m2 = 15, ma = 
4, m4 — 3, Ki — K2 — 0.1, Ti = 2, T2 = 1, T3 — 6,T4 — 7; (a) correlation plot between x and y laser 
outputs;(b)error xi — x dynamics; (c)error y\ — y dynamics. C is the correlation coefficient between 
X and y. Dimensionless units. 

FIG. 6. Numerical simulation of mutually coupled non-identical electro-optical laser systems, 
Eqs.(15-16): Partial generalized synchronization between x and y for Ki = K2 = 3.2(other pa- 
rameters as in Fig.5); (a) correlation plot between x and |/;(b)error xi— x dynamics;(c)error yi — y 
dynamics. C is the correlation coefficient between x and y. Dimensionless units. 
FIG. 7. Numerical simulation of mutually coupled non-identical electro-optical laser systems, 
Eqs.(15-16): Global generalized synchronization between x and y for Ki = K2 = 12(other pa- 
rameters as in Fig.5) (a) correlation plot between x and |/(b)error xi — x dynamics; (c) error yi — y 
dynamics. C is the correlation coefficient between x and y. Dimensionless units. 
FIG. 8. Numerical simulation of Eqs. (15-16), auxiliary systems' error dynamics yi — y, and xi — x 
for $0 = 7r/4,mi = 20, m2 — 15, ma = 4, m4 = 3, Ti = 2,T2 — 1,Ts — 6, r4 = 7.x laser system is 
more chaotic than y laser: The maximum conditional Lyapunov exponents (A) and ({>) 
against the coupling strength K for the auxiliary systems xi and x driven by y and yi and y 
systems driven by x, respectively. Dimensionless units. 
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FIG. 9. Numerical simulation of uni-directionally coupled non-identical electro-optical laser sys- 
tems, Eqs.(21-22): Retarded synchronization with the coupling delay time between x(solid line) 
and y (dotted hue) for ai = l,a2 = 18, $o = '''"/4, mi = = 7,m2 = m4 = 5,K = 17, Ti = 
3, T2 = 2, Ts = 5;C is the correlation coefficient between Xj-^ and y. Dimensionless units. 
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